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Abstract 

The cross sections for inelastic photoproduction of J /if) and tp' mesons have been meas- 
ured in ep colhsions with the ZEUS detector at HERA, using an integrated luminosity 
of 468 pb~^ collected in the period 1996-2007. The ip' to J /ijj cross section ratio was 
measured in the range 0.55 < z < 0.9 and 60 < < 190 GeV as a function of z 
and pt- Here W denotes the photon-proton centre-of-mass energy, z is the fraction of 
the incident photon energy carried by the meson and pT is the transverse momentum 
of the meson with respect to the beam axis. The J/tp cross sections were measured for 
0.1 < z < 0.9, QQ <W < 240 GeV and pr > ^ GeV. Theoretical predictions within the 
non-relativistic QCD framework including NLO colour-singlet and colour-octet contri- 
butions were compared to the data, as were predictions based on the /cp-factorisation 
approach. 
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1 Introduction 



The inelastic production of J/ip and of ip' has been studied for several years in hadron and 
electron-proton colliders and in fixed target experiments [1]. At HERA, the reactions 

ep eJ/i)X, (1) 

and 

ep ei^'X, (2) 

have been studied [2l[3] for low virtuality of the exchanged photon (photoproduction) in 
the range z < 0.9, where z denotes the fraction of the incident photon energy carried by 
the meson in the proton rest frame, thus excluding the diffractive process for which 2; ~ 1. 
In the HERA photoproduction regime, the production of inelastic J/ ip or ip' mesons arises 
mostly from direct and resolved photon interactions. In leading-order (LO) Quantum Chro- 
modynamics (QCD), the two processes can be distinguished; in direct-photon processes the 
photon enters directly into the hard interaction; in resolved-photon processes the photon 
acts as a source of partons, one of which participates in the hard interaction. The inelastic 
process in the photoproduction region is dominated by photon-gluon fusion. In this direct- 
photon process the photon emitted from the incoming electron interacts with a gluon from 
the proton to produce a pair of charm-anticharm quarks, cc, which then turn into the J / ip or 
the ip' mesons. When the cc pair emerges from the hard process with the quantum numbers 
of the mesons, the reaction is described in the framework of perturbative Quantum Chromo- 
dynamics (pQCD) by models such as the Colour Singlet (CS) model. In the Colour Octet 
(CO) model, the cc pair emerges from the hard process with quantum numbers different 
from those of the mesons and emits one or more soft gluons before turning into the physical 
meson state. Examples of direct-photon LO diagrams with a CS and a CO hard subprocess 
are shown in Fig. [TJ 

Full next-to-leading order (NLO) J /ip cross section predictions using only the direct-photon 
CS contributions have already been performed The non-relativistic QCD framework 

(NRQCD) [6] allows the evaluation of J /%p cross sections including direct and resolved photon 
processes with CS and CO contributions. The former contribution can be thought of as the 
first term of the NRQCD expansion and so it is an integral component of this theoretical 
formalism. Recently, the full computation was performed in the HERA photoproduction 
regime at the NLO level [7]. The numerical values of the CS and CO matrix elements 
were obtained from a global fit to hadroproduction, electroproduction and photoproduction 
inelastic J/ip data [7]. 
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J/ip cross sections have also been evaluated j8l|9] in the fcy-factorisation approach [10]. In 
this model, based on non-collinear parton dynamics governed by the CCFM evolution 
equations, effects of non-zero gluon transverse momentum are taken into account. Cross 
sections are then calculated as the convolution of unintegrated, transverse-momentum de- 
pendent gluon densities and LO off-shell matrix elements. Direct and resolved photon 
processes are included. The matrix elements are computed in the CS model. 

Measurements of the reactions ([T]) and ([2]) have been previously performed by the ZEUS 
collaboration |2], using an integrated luminosity of 38 pb~^, and by the HI collaboration j^, 
using an integrated luminosity of 165 pb"^. Total and differential cross sections were presen- 
ted as a function of various kinematical variables. The HI and ZEUS collaborations have 
also published a measurement of the J /ip helicity distribution [3l[I2], the ZEUS result was 
obtained using the full HERA luminosity. LO and NLO QCD predictions, as well as LO 
NRQCD calculations, were compared to the measurements. None of the calculations could 
describe the data in the whole kinematic range of the measurements. The data were shown 
to have the potential to reduce the large uncertainties in the phenomenological parameters 
used in the calculations. 

In this paper, measurements of reactions ([1]) and ([2]) are presented using a luminosity of 468 
pb~^. The J I'll) and ip' mesons were identified using the fi~^fi~ decay modes. 

The ip' to J/tp cross section ratio was measured in the range 60 < < 190 GeV and 
0.55 < z < 0.9 as a function of W, z and p^. Here W is the 7p centre-of-mass energy and 
•pT is the transverse momentum of the mesons with respect to the beam axis. The cross 
sections for inelastic J j'\p photoproduction as a function of p^, for different z ranges, and 
as a function of 2, for different "pj^ ranges, were measured in the range 60 < < 240 GeV, 
0.1 < 2 < 0.9 and Pt > 1 GeV. The momentum flow along and against the J/?/; direction of 
flight in the laboratory frame, as obtained from the charged tracks produced together with 
the Jji) in the range m <W < 240 GeV, 0.3 < ^ < 0.9 and 1 < pr < 10 GeV, was studied 
in order to shed further light on the production mechanisms. 



2 Experimental set-up 

The analysis presented here is based on data collected by the ZEUS detector at HERA 
in the period 1996-2007. In 1998-2007 (1996-1997), HERA provided electrorQ beams of 

^Here and in the following, the term "electron" denotes generically both the electron (e~) and the positron 
(e+). 
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energy Ee = 27.5 GeV and proton beams of energy Ep = 920 (820) GeV, resulting in a 
centre-of-mass energy of ^/s = 318 (300) GeV, giving an integrated luminosity of 430 (38) 

A detailed description of the ZEUS detector can be found elsewhere [TSlIT^ . A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) |15j . which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. Before the 2003- 
2007 running period, the ZEUS tracking system was upgraded with a silicon microvertex 
detector (MVD) [16]. In the following, the term "CTD-MVD track" denotes generically both 
the tracks measured in the CTD and (after 2002) in the CTD and MVD. 

The high-resolution uranium-scintillator calorimeter (CAL) |T^ consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeter^. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAG). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, as 
measured under test-beam conditions, were cr(E) / E = 0.18/-\/£' for electrons and (j{E) / E = 
0.35/\/^ for hadrons [E in GeV). The timing resolution of the CAL was better than 1 ns 
for energy deposits greater than 4.5 GeV. 

Muons were identified as tracks measured in the barrel and rear muon chambers (BMUON 
and RMUON) [TS]. The muon chambers were placed inside and outside the magnetised iron 
yoke surrounding the CAL. The barrel and rear inner muon chambers (BMUI and RMUI) 
covered the polar-angle regions 34° < 9 < 135° and 135° < 9 < 171°, respectively. 

The luminosity was measured using the Bethe-Heitler reaction ep — ?■ e7p with the lumin- 
osity detector which consisted of a lead-scintillator calorimeter jl9] and, after 2002, of an 
additional magnetic spectrometer [20] system. The fractional systematic uncertainty on the 
measured luminosity was 1.9%. 

^The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the proton- 
beam direction, referred to as the "forward direction", and the X axis pointing towards the centre of 
HERA. The coordinate origin is at the nominal interaction point. The polar angle, 9, is measured with 
respect to the proton-beam direction. The pseudorapidity is defined as ?7=-ln(tan |). 
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3 Event selection and kinematic variables 



The online and offline selections, as well as the reconstruction of the kinematic variables, 
closely follow the previous analysis |2]. 

Online, the BMUI and RMUI chambers were used to tag muons by matching segments in 
the muon chambers with CTD-MVD tracks, as well as with energy deposits in the GAL 
consistent with the passage of a minimum- ionising particle (m.i.p.). 

The different steps of the offline selection procedure are described in the following para- 
graphs. An event was accepted if it had two primary-vertex CTD-MVD tracks with invari- 
ant mass between 2-5 GeV. One track had to be identified in the inner muon chambers 
and matched to a m.i.p. cluster in the CAL. It was required to have a momentum greater 
than 1.8 GeV if it was in the rear region or a transverse momentum greater than 1.4 GeV 
if in the barrel region. The other track had to be matched to a m.i.p. cluster in the GAL 
and was required to have a transverse momentum greater than 0.9 GeV. Both tracks were 
restricted to the pseudorapidity region |?7| < 1.75. To reject cosmic rays, events in which 
the angle between the two muon tracks was larger than 174° were removed. 

In addition, events were required to have a calorimetric energy deposit larger than 1 GeV in a 
cone of 35° around the forward direction (excluding possible calorimeter deposits due to the 
decay muons). This requirement completely rejects exclusively produced J /il' mesons, ep — > 
epJ/ip. It also strongly suppresses the background from proton diffractive-dissociation, 
ep — eNJ/ip, because the low invariant mass hadronic system can often (but not always) 
escape along the outgoing proton direction without any activity in the FGAL. A reduction 
of the remaining background is achieved by requiring the events to have, in addition to the 
two decay muon tracks, at least one additional track with transverse momentum larger than 
250 MeV and pseudorapidity |?7| < 1.75. 

The ip' production in proton diffractive-dissociation processes with the decay chain J/ip^—^ 
/i"'"/i~) 7r+ TT^ was identified in the selected data sample. For the bulk of these events 
only four charged tracks are visible in the detector. For the events with a fi~^fi~ invariant 
mass, m^^, in the interval [2.85, 3.30] GeV, and with exactly two additional primary-vertex 
tracks of opposite charge, the total invariant mass m4 of the four tracks was evaluated. 
Events with a mass difference — m^^ within ±60 MeV of the nominal mass difference 
m^/ — mj/^ = 589 MeV [21] were discarded. This topology was tagged only in 1.2% of the 
overall selected J/ip sample and removed. 

These requirements effectively select inelastic J/i/) and ip' mesons. J/ip and ip' mesons from 
decays of b hadrons are also included in the data sample. 
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The kinematic region considered was defined by the inelasticity variable z and by the photon- 
proton centre-of-mass energy 

W'' = [P^q)\ (3) 

where P and q are the four-momenta of the incoming proton and the exchanged photon, 
respectively. It was calculated using 

W'' = 2E^{E-vz). (4) 

where [E — p^), the difference between the energy and the momentum along the Z axis, is 
summed over all final-state energy-flow objects |22] (EFOs) which combine the information 
from calorimetry and tracking. 



The inelasticity z = was determined as 



P-q 



(E-pz) ' 

where i/j can be either a J / i/j or a i/j' meson, p^ is the four-momentum of the tp and [E—pz)^ 
was calculated using the two tracks forming the ip. 

In order to reject deep inelastic scattering, events were required to have E — pz < 32 GeV. 
This restricts the virtuality of the exchanged photon, Q"^ = — g^, to ^ 1 GeV^, with 
a median of about 10~^ GeV^. The elimination of deep inelastic scattering events was 
independently confirmed by searching for scattered electrons in the CAL [23]; none was 
found. 

Table [T] summarises the various kinematic regions used for the presented measurements. 



4 Monte Carlo models 



The inelastic production of J/ip and ip' mesons was simulated using the HerwiG 6.100 |24] 
program, which generates direct photon events according to the LO diagrams of the photon- 
gluon fusion process, 7(7 — tpg. The processes are calculated in the framework of the CS 
model. The Herwig MC provides in general a good description of the data. To improve 
the agreement further, the p^ spectrum was reweighted to the data. The average weight 
of the MC events with pj- around 1 GeV is 0.85. The average weight for > 4 GeV is 
instead 1.8. 
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Diffractive production of J/iIj and ip' mesons with proton dissociation was simulated with 
the Epsoft j25] MC generator, which was tuned to describe such processes at HERA |26) . 

The Pythia 6.220 MC generator [27] was used to generate J/ip and Xc states from the 
resolved-photon process, with LO matrix elements computed in the CS model. The gen- 
erator cross sections for the J/ip and Xc states are very similar. For the generation of the 
Xci(l-P) and Xc2(l-P) mesons, only the J /ip decay channel was considered. The final state 
photon is at low energy, 0(400) MeV, basically indistinguishable from the remaining had- 
ronic activity of the event. Hence the effective resolved-photon J/il) contribution can be 
thought of as due to the genuine resolved-photon component plus the Xc feed-down. The 
resolved ip' contribution was neglected due to the small resolved-to-direct cross section ratio 
and to the additional reduction due to the ip' — > J/ipX branching ratio. 

The Pythia MC was also used to generate the production of J / ip and ip' mesons originating 
from h hadron decays, mostly from 5-mesons. The following beauty-quark production 
processes were generated (according to the Pythia notation): direct, resolved, 7 and proton 
excitation. The beauty-quark mass was set to 4.75 GeV and the branching ratios of the h 
hadrons to J / tp and tp' were set to the corresponding PDG [21] values. 

All generated events were passed through a full simulation of the ZEUS detector based on 
Geant 3 [28]. They were then subjected to the same trigger requirements and processed 
by the same reconstruction program as the data. 

5 Signal determination and cross sections calculation 

The invariant-mass spectrum of the muon pairs measured in the phase space region used in 
the determination of the ip' to J/ip cross section ratio, QO < W < 190 GeV and 0.55 < z < 
0.9, is shown in Fig. [2l A non-resonant background contribution, mostly due to hadrons 
misidentified as muons, is also visible. This contribution was estimated by fitting the product 
of a second-order polynomial and an exponential function to the region 2-2.75 and 3.8-5 
GeV, outside the J/ip and ip' invariant-mass window. The number of J/ip events was 
obtained by subtracting the number of background events, estimated from the fit procedure, 
from the total number of events inside the J/tp invariant-mass window, 2.85-3.3 GeV. This 
procedure resulted in 11295 ± 114 J/tp events. The same procedure applied to the ip' 
invariant-mass window, 3.55-3.8 GeV, gave 448 ± 34 events. 

Applying the same procedure to the phase space region used for the differential J/ip cross 
section measurements, 60 < < 190 GeV, 0.1 < 2; < 0.9 and pt > I GeV, 12671 ± 
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161 J/il) events were found. The fitting procedure described above was performed for each 
measurement bin presented in this paper. 

The cross section for any observable, O, was computed for each bin, i, using correction 
factors, Ci{0), defined as Ci{0) = Nf''{0)/Nl''\0), where Nf'^O) is the number of 
events generated with the Herwig MC and Nl'^^iO) is the number of the events recon- 
structed by the standard analysis chain. The factors Ci{0) take into account the overall 
acceptance including the geometrical acceptance and the detector, trigger and reconstruction 
efficiencies. They also take into account bin-to-bin migrations. 

For 0.9 < z < 1, the events are largely diffractive. Therefore, the analysis of inelastic J /il) 
production was restricted to the region 0.1 < z < 0.9. In order to further suppress diffractive 
events, the transverse momentum of the J/ip mesons had to fulfill > 1 GeV. The 
remaining contamination was estimated by fitting the relative fractions of non- diffractive and 
diffractive events to the data ^-distribution, using the HerwiG and EPSOFT MC simulations 
as templates. From this fit, the overall diffractive background contribution for 0.1 < z < 0.9 
is 4.6 ± 1.6%. 

In Fig. |3] the Herwig and EPSOFT MC mixture, in the kinematic region < W < 240 
GeV, 0.3 < z < 0.9 and p^p > 1 GeV, is compared to the data: a reasonable description 
is found. The region 0.1 < z < 0.3 was removed because not diffractive background is 
present at low z. The estimated diffractive background was subtracted bin by bin from the 
measured differential cross sections. 

The cross sections measured in this analysis include also contributions from resolved-photon 
processes and from decays of beauty hadrons. Inelastic J/ip production via the resolved- 
photon process has not been measured explicitly up to now in the photopro duct ion regime. 
QCD predictions, as well as the Pythia MC simulation described in Section |U indicate 
that this contribution is largest at low z values. For z < 0.1, the expected size of this 
contribution can be larger than the direct-photon component. However, for z > 0.1, the 
resolved-to-direct photon production ratio is expected to be small. Since the acceptances 
obtained from the Herwig and Pythia MC simulations are similar, the Herwig MC alone 
was used for the overall acceptance corrections. 

The contribution to the measured cross sections due to J/ ip originating from B meson decays 
was estimated using the inclusive beauty Pythia MC sample described in Section |H The 
simulation predictions were scaled by a factor 1.11 according to the recent ZEUS measure- 
ment [29] of beauty photoproductioij^l. This leads to the estimation that on average 1.6% 

^The scaling factors obtained in the measurements [5511301 vary between 1.11 and 1.84 
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of the observed J/il) mesons originated from beauty hadron decays. The largest relative 
contribution, 4.5%, is in the kinematic region 0.1 < z < 0.3 and 1 < < 2 GeVl This 
component is not subtracted from the measured cross sections. 



6 Systematic uncertainties 

For all the measured quantities, the following sources of systematic uncertainties were in- 
vestigated (their effects on the measured cross sections are given in parentheses): 

• muon trigger and reconstruction efficiencies: the BMUI and RMUI muon chamber 
efficiencies were extracted from the data using muon pairs from elastic J/ip events 
and from the process 77 — fi~^fi~. These efficiencies take into account the full muon 
acquisition chain, from the online to the offline level and are known with a ±5% 
uncertainty (5% uniformly distributed in pj- and z); 

• hadronic energy resolution: the W and z resolutions are dominated by the hadronic 
energy resolution affecting the quantity {E — pz). The hadronic {E — pz) resolution 
in the MC was smeared event by event by ±20%, a conservative upper limit of a pos- 
sible systematic difference between data and MC. This gave only small cross sections 
variations (< 5%); 

• HerwiG MC Pt spectrum: the pt spectrum of the J ftp mesons in the Herwig MC 
simulation was varied within ranges allowed by the comparison between data and 
simulation and the correction factors were re-evaluated (< 5%); 

• J/ip helicity distribution: the J/ip helicity distribution can be described by two para- 
meters A and u [31]. In the Herwig MC these are set to zero. According to the direct 
measurement of the helicity parameters performed by ZEUS [T^, all data points lie 
within the region of the X-u plane defined by |A| < 0.5 and \u\ < 0.5 with only a 
mild Pt or z dependence. Hence, as a systematic check, the Herwig MC was re- 
weighted varying independently A and u in the range ±0.5 and the correction factors 
were re-evaluated (5 — 10% depending on the px and z region); 

• diffractive simulation: the EPSOFT MC simulation parameters were varied within 
ranges allowed by the comparison between data and the EPSOFT MC simulation in 
the region 0.9 < z < 1. The diffractive background was re-evaluated (< 5% at high z 
and low pt, negligible elsewhere); 
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• diffractive subtraction: the relative fraction of inelastic and diffractive processes, as 
represented by the Herwig and EPSOFT MC, was fixed by the procedure described 
in Section [5l It is known to a precision limited by the number of J/ip events in the 
data and the process modeling by the MCs. The relative fractions were varied within 
ranges allowed by the comparison between data and simulation (up to 10% at high z 
and low pti negligible elsewhere); 

• invariant-mass window: the m^+^- invariant-mass window used to estimate the num- 
ber of J/ip events above the non-resonant background was enlarged to [2.8,3.35] GeV 
and tightened to [2.9,3.3] GeV. For the ip' to J/il) cross section ratios, similar mass 
window variations were also applied for the ip' signal (generally < 5%, up to 10% 
at low z values where the number of expected and observed events is small and the 
non- resonant background is largest); 

• additional track cut: the requirement of three tracks, including the two J/ if) decay 
muons, with transverse momentum larger than 250 MeV and pseudorapidity \ri\ < 
1.75, was replaced by the requirement of five tracks with transverse momentum larger 
than 125 MeV, in the same pseudorapidity range. With this stronger requirement the 
diffractive J/ip background and the diffractive ip' contribution via the cascade decay 
J/ip{^ fj''^fi~) TT^ IT" are expected to vanish. Furthermore, a change in the overall 
multiplicity cut allows a test of how well the MC model reproduces the data in this 
respect. The MC mixture gives a fair description of the track multiplicity observed 
in the data. The cross sections were re-evaluated with the harder multiplicity cut 
(generally < 5%, up to 20% in some bins at low z and high z high px)- 

All of the above individual sources of systematic uncertainty were added in quadrature. 
The following sources would result in an overall small shift of the cross sections: 

• the integrated luminosity determination gave an uncertainty of ±1.9%; 

• the J/ifj — )■ fi~^fi~ branching ratio, 5.93 ± 0.06% gave an uncertainty of ±1%. 
They were not included. 



9 



7 Results 



7.1 ip' to J/ip cross section ratio 

The Ip' to J/ip cross section ratio was measured using the rates of ip' — t- fi^fi~ and J/ip ^ 
/i+yU". The ratio was determined in the region 60 < < 190 GeV, 0.55 < z < 0.9. The 
Pt > I GeV requirement was removed to maximise the available statistics. An increase 
of the diffractive background is expected. But under the assumption that this background 
contribution will be the same for the ip' and J/ip mesons it will cancel in the cross section 
ratio. The range 190 < W < 240 GeV and 0.1 < z < 0.55 was not included because the ip' 
peak was not visible in this high W and low z region. The ip' to J/ ip cross section ratio was 
computed in bins of W, z and px from 



where, for the considered bin i, N^^ {Nf'^) is the number of J/ip {ip') events observed, C^^ 
(Cf^) is the correction factor (see Section [5]) computed using the Herwig MC, Br^ {Br^') 
is the J/ip (ip') muonic branching ratio and Br' is the ip' ^ J/ip X branching ratio. The 
values used are Br^' = 5.93%, Br>'' = 0.77% and Br' = 59.5% |2I]. With this technique, 
the cross section ratio was corrected for the ip' ^ J/ip (— ^~^f^~) X cascade decay. 

Since NLO predictions are not available for ip', only the LO CS model expectations can be 
compared to the data. In the CS model, the underlying production mechanism is the same 
for J/ip and ip', hence all cross section ratios should be largely independent of the kinematic 
variables. Using the values of Br^ and Br^' given above, the expected ratio is 0.25 jl]. 
Since the NLO corrections, though being large, should be similar for J/ip and ip', the ratio 
at NLO is not expected to differ significantly from that at LO. 

The results, shown in Fig. [Hand listed in Tabled are dominated by the statistical uncertain- 
ties while most of the systematic uncertainties cancel in the ratio. The LO CS predictions 
agree reasonably well with the data. 

7.2 J/i/j differential cross sections 

The J/ip differential cross sections presented here include the inelastic ip' feed-down via 
the decay ip' — )■ J/ip (— )■ fi^fJ^~) X and the contribution from b hadron decays. The ip' 
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feed-down contributes about 15% and the b hadron decays 1.6% (see Section [5]). The W 
range of the differential cross sections is 60 < < 240 GeV. 

The differential cross sections da/dp^ were measured in the range 1 < p\ < 100 GeV^ 
for different z ranges. The results are listed in Table [3] and shown in Figs. [S] and O The 
predictions of a NRQCD calculation are compared to the data in Fig. [5] and those based 
on the /cT^-factorization approach [H] in Fig. ElH 

The differential cross sections da /dz were measured in the range 0.1 < z < 0.9 for different 
Pt ranges. The results are shown in Figs. [7] and [S] and listed in Table HI 

The present measurements are in agreement with the results obtained by the HI collabor- 
ation [3] except in the region z > 0.6 and px > 3 GeV where the ZEUS cross sections are 
above the HI measurements. 



7.2.1 Comparison of NRQCD calculation 



In Fig. |5]a prediction [7] performed in the NRQCD framework including direct and resolved 
photon processes is compared to the measured da/dp^. The hard subprocesses take into 
account both CS and CO terms to NLO. The square of the renormalisation and factorisation 
scales used is 4 ■ + p^, the charm quark mass, m^, is set to 1.5 GeV and the strong 
coupling constant, ^^(M^), to 0.118. The NRQCD scale, connected to the colour-octet 
terms, is set to m^. The CS contribution alone predicts cross sections significantly below 
the datcfl and fails to describe the data in all z regions shown here. Including CO terms 
give a dramatic improvement and leads to a rough agreement with the data. In general the 
calculation reproduces the steep drop of da/dp"^ with p^, however, in the intermediate z 
range, 0.3 < z < 0.75, the prediction rises less steeply than the data towards the smallest 
values of p"^. 

In Fig. [7] the NRQCD predictions described above are compared to the measured da/dz. 
The predictions rise too steeply with z compared to the data, for all the p^ ranges. 

*Both the NRQCD and the fcr-factorisation calculations do not include tp' feed-down and b hadron decays, 
however these expected contributions are small compared to the uncertainties of the calculations. 

^ The NLO CS predictions [4] shown in the previous publication |2] were the first performed and used 
extreme values for the renormalisation and factorisation scales, with the effect of artificially increasing 
the normalisation of the predicted cross sections |32| . 
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7.2.2 Comparison of A; j^— factorisation approach 

In Fig. |6] a prediction [9] performed in the fc-r-factorisation approach is compared to the 
measured da/dp^. The matrix elements are computed in the CS model using rric = 1.5 GeV 
and as{Mz) = 0.1232. In the numerical calculation, the renormalisation and factorisation 
scales squared are set to tu'^j/^ +Pt ^^'^ ^ + Qt^ respectively, where s is the four-momentum 
squared of the hard subprocess and Qt is the transverse momentum of the initial parton. 
The unintegrated CCFM parton density |33| was selected. Using different sets of parton 
densities leads to changes in the prediction that are small with respect to the effects of 
scale variations already shown in Fig. O Thus this source of theoretical uncertainties was 
neglected. The /cr^factorisation prediction, with the values of rric and as given above, 
provides a better description of the data than the NRQCD model. 

The above /cT^factorisation predictions are compared to the differential cross sections da/dz 
in Fig. [SI Here too the description is better than that of the NRQCD model. Note however 
that the /cr^factorisation model prediction suffers from large theoretical uncertainties, in 
particular at low p^- 

7.3 Momentum flow along and against the J/i/j direction 

As pointed out by Brambilla et al. pOj, the different colour flow in CS and CO hard subpro- 
cesses is expected to translate into different properties of the hadronic final state. In the 
photoproduction regime, the transverse momentum of the incoming photon is negligible. 
Thus in the CS model (see Fig. [T](a)), at LO the J/tp and the final state gluon are expected 
to be back to back. Hence, in this model, the momentum flow along the J/ip direction, P^iong , 
is expected to be small. The momentum flow against the J/ ip direction, Pagainst ? should in- 
stead be driven by the hadronisation of the gluon. In the CO framework (see Fig. [U (b)), 
no substantial difference is expected for Pagainst , compared to the CS framework. Instead, a 
contribution due to the soft gluons emitted by the cc pair forming the physical J/ip state 
should be present. Hence, Pagainst is again sensitive to gluon fragmentation while P^iong can 
shed light on the CO dynamics. As NRQCD framework MC generators are not presently 
available for ep collisions, only predictions of the CS model Herwig MC are compared to 
the data. 

The momentum flow analysis was performed for different px ranges. All track quantities 
described in the following were measured in the laboratory frame at the reconstruction level. 
Only primary vertex tracks with px > 150 MeV and \ri\ < 1.75 were selected. The J/ip 
decay muon tracks were discarded. For each track whose component of momentum along 
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the J/il) direction in the laboratory frame was positive, the component was included in 
Paiong • If it was negative, it was included, in absolute value, in Pagainst • The data were 
restricted to z > 0.3 where the signal to background ratio is highest. The W and ranges 
were < W < 240 GeV and 1 < < 10 GeV, respectively. The residual non-resonant 
background was subtracted for both Pagainst and Paiong variables using the shapes measured 
in the J/ip side bands region and the normalisation obtained from the signal extraction 
procedure described in Section O 

The Pagainst (Paiong ) distrlbutiou, normalized to one, is shown in Fig. [H] f|TU|) . The prediction 
obtained from the Herwig MC simulation (including detector simulation) is also shown. 
The Pagainst distrlbutiou of the MC simulation shows a softer drop from the first to the second 
momentum bin than that of the data. This situation is reversed for the higher momenta 
values where Herwig predicts a steeper decrease than that observed in the data. This 
behavior is seen for all Pt regions. 

For the Paiong distribution, shown in Fig. [101 a better agreement is found between the Herwig 
MC prediction and the data. 

8 Conclusions 

A measurement of the inelastic photopro duct ion of J / and ip' mesons at HERA was presen- 
ted. The tp' to J ftp cross section ratio was measured as a function of several kinematical 
observables. The constant value of 0.25 predicted by the LO CS model is in reasonable 
agreement with the data. 

Double differential cross sections of inelastic J/ip photoproduction were measured. A LO 
calculation [H] using CS terms alone gives, within large normalisation uncertainties, a 
good description of the differential cross sections. However, for a better comparison with 
the data, a reduction of the theoretical uncertainties is very important. 

A recent NLO calculation [7], using CS and CO terms in the collinear approximation, gives 
a rough description of the double differential cross sections. The same calculation with only 
CS terms is in strong disagreement with the data. This leads to the conclusion that CO 
terms are an essential ingredient for this particular model. 

Predictions of the Herwig MC, which includes only CS processes, were compared to the 
measured momentum flow along and against the J/ip direction. Herwig reproduces the 
fall off of the momentum distribution against the J/ ip direction as the momentum increases 
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but fails to describe the exact shape of this distribution. A better description is obtained 
along the J / direction. 
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ip to J/ip cross section ratio: kinematic range 


60 <W < 190 GeV 


PT>0 GeV 


0.55 < z <0.9 


Differential cross sections: kinematic range 


60 <W < 240 GeV 


Pt > 1 GeV 


0.1 < z < 0.9 


Momentum flow: kinematic range 


60<W < 240 GeV 


1< PT < 10 GeV 


0.3 < z < 0.9 



Table 1: The different kinematic regions used in the measurement of the ip' to J/ip cross 
section ratio, J/ip differential cross sections and momentum flow along and against the J/ip 
direction. 



Pt range 
(GeV) 


{pt) 
(GeV) 


a{i'')/a{J/i') 


0.0-1.0 
1.0-1.75 
1.75-5.0 


0.63 
1.35 
2.68 


0.262 ± OMStoml 
0.317 ±0.049t^:^J^ 
0.263 ±0.04li°:°^° 


W range 
(GeV) 


(W) 
(GeV) 


a{i'')/aiJ/i') 


60-95 
95 - 120 
120 - 190 


81.42 
108.03 
149.11 


0.368 ± 0.054j:°;^^2 

0.409 ± 0.057^^15 
0.218 ±0.040t|]:o?^ 


z range 


(^> 


aii'')/a{J/i') 


0.55-0.70 
0.70-0.80 
0.80-0.90 


0.62 
0.75 
0.85 


0.250 ± 0.043t^;J^i5 
0.289 ± OMOtoml 
0.344 ± 0.054lg;^^^ 



Table 2: Cross section ratio of ip' to J/ip as a function of pt, W and z in the kinematic 
region 60 < W < 190 GeV and 0.55 < z < 0.9. In the quoted ratios, the first uncertainty is 
statistical and the second is systematic. 
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z range 


p|, range 




da / dpip 


da{b — s- J l^)ld%>\ 




(GeV2) 


(GeV2) 


(nb/GeV2) 


(nb/GeV2) 


0.10 - 0.30 


1.0 - 2.0 


1.46 


1.03±0.13to;[5 


0.05 




2.0 - 3.0 


2.47 


0.86±0.12l°-:t? 


0.04 




3.0 - 4.5 


3.67 


0.410 ± 0.079lo o5| 


0.029 




4.5 - 7.0 


5.64 


0.127 ±0.047lo o2? 


0.018 




7.0 - 10.0 


8.37 


0.052 ± O.OSOto oog 


0.012 




10.0 - 14.0 


11.62 


0.056 ± ^-OntoMd 


0.007 




14.0 - 20.0 


16.34 


0.0329 ± 0.008lto oogg 


0.0035 




20.0 - 40.0 


26.50 


0.0069 ± 0.0018j:o ooJg 


0.0012 




40.0 - 100.0 


56.69 


0.00092 ± 0.00037to;ooo26 


0.00013 


0.30 - 0.45 


1.0 - 2.0 


1.47 


1.32±0.10+oig 


0.02 




2.0 - 3.0 


2.45 


0.823 ± 0.08ll°;?|5i 


0.018 




3.0 - 4.5 


3.70 


0.492 ± 0.060lo o75 


0.013 




4.5 - 7.0 


5.64 


0.190 ±0.032i{j;!j2| 


0.010 




7.0 - 10.0 


8.35 


0.111 ±0.019^0.0" 


0.006 




10.0 - 14.0 


11.77 


0.062 ± O.Olllo gJ? 


0.004 




14.0 - 20.0 


16.49 


0.0349 ± 0.0052l[5;°°3° 


0.0021 




20.0 - 40.0 


27.96 


0.0065 ± Qmi2tl °QQQfi 


0.0007 




40.0 - 100.0 


54.05 


0.00095 ± 0.00019l||!;5]g]^j7 


0.00009 


0.45 - 0.60 


1.0 - 2.0 


1.45 


2.20 ± 0.09^0 25 


_ 




2.0 - 3.0 


2.47 


1.38±0.08l°;i5 


- 




3.0 - 4.5 


3.69 


0.84 ± O.OSio io 


- 




4.5 - 7.0 


5.65 


0.424 ± 0.029lo ggg 


- 




7.0 - 10.0 


8.35 


0.249 ± 0.017to.o29 


- 




10.0 - 14.0 


11.79 


0.121 ± 0.010lo oi3 


- 




14.0 - 20.0 


16.60 


0.0505 ± 0.0048lE!;5j^*^ 


0.0007 




20.0 - 40.0 


26.70 


0.0106 ± 0.0011155;°°°^ 


0.0004 




40.0 - 100.0 


55.86 


0.00122 ± 0.00020to;oooi5 


0.00003 


0.60 - 0.75 


1.0 - 2.0 


1.45 


2.80±0.10l°;32 


- 




2.0 - 3.0 


2.47 


2.07±0.09lo23 


- 




3.0 - 4.5 


3.70 


l.lOitO.OSlo 


- 




4.5 - 7.0 


5.60 


0.680 ± 0.030lo og4 


- 




7.0 - 10.0 


8.38 


0.286 ± 0.017tg o3g 


- 




10.0 - 14.0 


11.93 


0.153 ± O.OlOio ols 


_ 




14.0 - 20.0 


16.92 


0.0532 ± 0.0044lE5;{jogJ 


- 




20.0 - 40.0 


27.0 


0.0123 ±0.001llo;o°" 


- 




40.0 - 100.0 


55.77 


0.00112 ± 0.000181o;qoq22 


- 


0.75 - 0.90 


1.0 - 2.0 


1.45 


2.39±0.13to35 


- 




2.0 - 3.0 


2.45 


1.77±0.1llO-«3 


- 




3.0 - 4.5 


3.66 


1.17±0.07l5j;i^ 


- 




4.5 - 7.0 


5.64 


0.716 ± 0.039lo o|2 






7.0 - 10.0 


8.31 


0.369 ± 0.023to o45 






10.0 - 14.0 


11.77 


0.166 ± O.012l:5j;[',2o 






14.0 - 20.0 


16.66 


0.0650 ± 0.0058to;oogf 






20.0 - 40.0 


26.22 


0.0139 ± 0.0013to o°25 






40.0 - 100.0 


54.0 


0.00093 ± 0.00018tg;°oo?g 





Table 3: Measured J differential photoproduction cross sections in the kinematic region 
0.1 < z < 0.9 and QO < W < 240 GeV as a function of the squared transverse momentum of 
the J/ip mesons in bins of inelasticity z. IriJ^e quoted cross sections, the first uncertainty 
is statistical and the second is systematic. The bin center values {p^) and the expected, but 
not subtracted, beauty contribution (estimated through the Pythia MC) are also given in 
the table. The beauty contribution is only given when its value is above 1% with respect to 
the corresponding measured differential photoproduction cross section. 



Pt range 
(GeV) 


z range 




da/dz 
(nb) 


da{b — > J/'tJj)/dz 
(nb) 


1.0-2.0 


0.10-0.30 
0.30-0.45 
0.45 - 0.60 
0.60-0.75 
0.75 - 0.90 


0.21 
0.37 
0.52 
0.67 
0.82 


11.5± 1.0l}j 
17.3± l.Ot^J 

29.9 ± 0.9tli 

40.2 ± 1.0^4 5 
36.6± 1.2t'^| 


0.6 
0.4 


2.0-3.0 


0.10-0.30 
0.30-0.45 
0.45 - 0.60 
0.60-0.75 
0.75 - 0.90 


0.21 
0.37 
0.52 
0.67 
0.82 


1.94±0.78t^:!^5 
6.42±0.7lt°:^^ 
12.4 ± 0.61:} j 

18.6±0.6l2:2 
19.4±0.8;2:3 


0.43 
0.28 
0.8 


3.0-4.5 


0.10-0.30 
0.30-0.45 
0.45 - 0.60 
0.60-0.75 
0.75-0.90 


0.20 
0.38 
0.52 
0.68 
0.82 


2.55 ±0.47+^:^^ 

6.61 ± 0.37l|]:?^ 
7.79 ± 0.35t!}:^5 
9.29 ± 0.46t'j':^5 


0.30 
0.26 
0.08 


> 4.5 


0.10-0.30 
0.30-0.45 
0.45 - 0.60 
0.60-0.75 
0.75 - 0.90 


0.21 
0.38 
0.52 
0.67 
0.82 


1.01 ± 0.20t^:i|^ 
1.31±0.18t|]:}^ 
1.98±0.17to:i6 

2.ii±o.i6tg:i^ 

2.16±0.18t°:^^ 


0.16 
0.24 
0.08 



Table 4: Measured J ftp differential photoproduction cross sections in the kinematic region 
Pt > I GeV and 60 < < 240 GeV as a function of the inelasticity z in bins of transverse 
momentum of the J / ip meson. In the quoted cross sections, the first uncertainty is statistical 
and the second is systematic. The bin center values (z) and the expected, but not subtracted, 
beauty contribution (estimated through the Pythia MG) are also given in the table. For 
further details see Tablel^ 
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Figure 1: Examples of direct photon-processes at leading- order in fa) the colour- singlet 
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Figure 2: Invariant-mass distribution, rrinn, in the kinematic region 0.55 < 2; < 0.9 and 
60 < W < 190 GeK The continuous line shows the estimated background contribution (for 
further details see the text). The right insert highlights the ip' mass peak. 
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Figure 3: J/tp events fraction measured in the kinematic region 0.3 < z < 0.9, 60 < < 
240 GeV and > 1 GeV as a function of (a) the polar angle 6f^ of the muon tracks, (b) 
W , (c) the inelasticity z and (d) the J/ip px- The data are shown as points. The error 
bars are the statistical uncertainties. The sum of the Herwig and EPSOFT MC predictions, 
according to the relative fraction described in the text and normalised to the data are also 
shown (continuous lines). The EPSOFT MC component is shown separately (dashed lines). 
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Figure 4: xp' to J /if) photoproduction cross section ratio measured in the kinematic region 
0.55 < z < 0.9 and 60 <W < 190 GeV as a function of (a) W, (h) the inelasticity z and (c) 
Pt- The data are shown as points. The inner error bars are the statistical uncertainties, while 
the outer error bars show the statistical and systematic uncertainties added in quadrature. 
The leading-order colour-singlet model expectation (horizontal lines) is also shown. 
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Figure 5: Differential cross sections da/dp^ measured in 5 different z ranges. The meas- 
urement is performed in the kinematic region 60 < TV < 240 GeV and > 1 GeV. The 
data are shown as points. The inner (outer) error bars represent the statistical (total) un- 
certainties. The solid lines show the NLO CS+CO (BK) prediction JTj obtained in the non- 
relativistic QCD framework. The uncertainties are indicated by the band. The colour-singlet 
model contribution is presented separately as the dashed lines. 
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Figure 6: Differential cross sections da/dp^ measured in 5 different z ranges. The meas- 
urement is performed in the kinematic region 60 < < 240 GeV and px > 1 GeV. The 
data are shown as points. The inner (outer) error bars represent the statistical (total) uncer- 
tainties. The solid lines show the kx -factorisation (BLZ) prediction JB^. The uncertainties 
are indicated by the band. 
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Figure 7: Differential J/ip cross sections da/dz measured in 4 different px ranges. The 
measurement is performed in the kinematic region 60 < W < 240 GeV and 0.1 < z < 0.9. 
The data are shown as points. The inner (outer) error bars represent the statistical (total) 
uncertainties. The solid lines show the NLO CS+CO (BK) prediction obtained in the 
non-relativistic QCD framework. The uncertainties are indicated by the band. 
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Figure 8: Differential J/ip cross sections da/dz measured in 4 different px ranges. The 
measurement is performed in the kinematic region 60 < W < 240 GeV and 0.1 < 2; < 
0.9. The data are shown as points. The inner (outer) error bars represent the statistical 
(total) uncertainties. The solid lines show the kx -factorisation (BLZ) prediction JB^. The 
uncertainties are indicated by the band. 
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Figure 9: Momentum flow against the J/ip direction of flight in the laboratory frame, 
Pagainat , for different px ranges. The distributions are normalized to unity and are not 
corrected for detector acceptance. The measurement is performed in the kinematic region 
60 < W < 240 GeV and 0.3 < z < 0.9. The data are shown as points with error bars in- 
dicating their uncertainties. The predictions obtained from the Herwig MC are also shown 
as rectangular shaded boxes. The height of these boxes represents the uncertainties of the 
prediction. 
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Figure 10: Momentum flow along the J/ip direction of flight in the laboratory frame, Paio-ng , 
for different pj- ranges. The distributions are normalized to unity and are not corrected for 
detector acceptance. The measurement is performed in the kinematic region 60 < W < 240 
GeV and 0.3 < z < 0.9. The data are shown as points with error bars indicating their 
uncertainties. The predictions obtained from the Herwig MC are also shown as rectangular 
shaded boxes. The height of these boxes represents the uncertainties of the prediction. 
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